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ABSTRACT 

Recently, I have defined the so called PDF's (prime distribution factors) which 
govern the distribution of prime numbers of the type p,p + ai being all primes up to 
some number n. It was shown that the PDF's are expressible in terms of the basic 
PDF's which are defined as Oj — dj or a{ being composed of primes which are less or 
equal to the number of primes. For example, p,p + 2 (twin primes), or p,p + 2,p + 6 
being all primes (basic triplets). We give here a conjecture for the number of 
basic prime PDF's in terms of Hardy-Littlewood numbers, thus completing the 
determination of PDF's. These conjectures are supported by extensive calculations. 



The distribution of primes has been fascinating since the work of Gauss and 
Riemann. In ref. [1] the following observation was made. Define by N(ai,n) the 
number of sequences of p, p + ai, i = 1 ... m — 1 being all primes, up to the integer n. 
Then the following limit exists and C(a,i) is termed the "PDF" (prime distribution 
factor) : 

n 
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for large enough n. The PDF's were shown to be expressible in terms of the basic 
PDF's by a product with some rational number. For example, the PDF for p and 
p + n (n even) is given by 

C(p,p + n) = C(p,p + 2) J] i^l. (2) 

q>2 prime 
q\n 



The basic PDF's are defined by the requirement that aj — aj and d{ are all 
products of primes less or equal to m. For example: p, p+2, for twins, p, p + 2,p+6, 
for triplets, p, p + 2,p + 6,p + 8 for quadruplets, p, p + 2, p + 6,p + 8,p + 12 for 
quintuplets, add p + 18 for sextuplets, etc. 

A conjecture for prime twins, p,p + 2 is known before, see, e.g., [2]. It is given 

by 

C(2) = 2c 2 , (3) 



where 



p>2 
p prime 



(P-1) 

■rime 

and is numerically C(2) = 1.32032... 
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Our generalization to this is to use Hardy-Littlewood numbers, defined by 



Cm ~ 11 (p-l)m ' W 

p prime 

The first few HL (Hardy-Littlewood) numbers are: C2 = 0.6601618, C3 = 0. 6351663. ...| 
C4 = 0.3074948..., C5 = 0.409874..., etc. See, e.g., ref. [3], where these numbers are 
calculated, and references therein. Our main observation is that the basic PDF's 
are given by the HL numbers times a "simple" rational number: 

Conjecture: The basic PDF's are given by HL numbers times a rational num- 
ber. This number is 2 for twins, 9/2 for triplets, 27/2 for quadruplets. 

Thus, we have for the basic PDF's the formula: 



C(2) = 2c 2 = 1.32032... (6) 

C(3) = 9c 3 /2 = 2.858248... (7) 

C(4) = 27c 4 /2 = 4.1511808... (8) 

To check these conjectures we made an extensive calculation of the number of 
basic m-plets. The program we used was written in MATLAB and is listed in the 
appendix. 

For twin primes we find up to 10 10 , 27, 412, 673 twins, p and p + 2 being primes. 
This gives the PDF 1.32038... which is 4.5 x 10~ 5 off the limit C(2), eq. (6). For 
up to n = 2 x 10 10 we find 4, 942, 554 basic triplets, of the form p, p + 2,p + 6, 
giving the PDF 2.85768 which is 2 x 10~ 4 off C(3), eq. (7). For basic quadruplets, 
p,p+2,p+6,p+8, we find 898, 998 up to the number 7x 10 10 , which gives the PDF, 
4.1503.., which is 2 x 10~ 4 off C(4), eq. (8). This supports the above conjecture. 
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For basic quintuplets we were unable to determine the rational number, which 
seems not as simple as for m < 4. For n = 4 x 10 11 we find 370, 502 basic 
quintuplets, which gives the PDF, 10.1193 or C(5) ~ 24.6888 x C5. Thus, we know 
that the rational factor is 24.6888 approximately, but it is not enough to form a 
conjecture about it. 

To summarize, we have a conjecture about the number of all prime multiplets, 
expressing them in terms of Hardy-Littlewood numbers, times some known rational 
coefficients. These completes the calculation of the number of prime multiplets up 
to quadruplets. For higher multiplets, we were unable to determine the coefficients, 
which seems not so simple, and it is left to further work. 
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APPENDIX 

Attached are two MATLAB functions used to calculate the number of basic 
prime multiplets. We used the function, "primel2", to determine the number of 
primes between nl and n2, employing the sieve method. Further below is a sample 
function, "pairs5", which calculates the number of basic quintuplets up to the 
integer STEPx mm, which works by looking on consequent primes from a list 
generated by "primel2", and determining their differences. 

function p = primel2(nl,n2) 

%PRIMEL2 Generate list of prime numbers. 

% PRIMEL2(N1,N2) is a row vector of the prime numbers less than or 
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% equal to N2 and bigger or equal to Nl. 

pl=primes(sqrt(n2)); 

lpl=length(pl); 

qq=nl:n2; 

ifnl==l; qq(l)=0; 

end; 

for mm=l:lpl 
qw=pl(mm); 

sw2=nl-mod(nl-l,qw)+qw-l; 
if sw2==qw ;sw2=sw2+qw; 
end; 

qq(sw2-nl+l:qw:n2-nl+l)=0; 
end; 

p=qq(qq > 0); 

% 
% 

% Below is the sample function, "pairs5", which calculates the number 
% of basic quintuplets, up to the integer mm x step. The reason of doing 
% this in "batches" is to prevent exceeding the available memory. 
% We typically take step=10 7 . 

% 

function np = pairs5(mm,step) 
np=0; 

for jj=l:mm; 
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nl=step*(jj-l)+l; 

n2=step*jj; 

if nl > 1; nl=nl-ll; 

end; 

pp=primel2(nl,n2); 

np=np+sum([0 pp 0]-[0 pp]==2 & [0 pp 0]-[0 pp 0]==4 ... 

& [0 pp 0]-[0 pp 0]==2 & [pp 0]-[0 pp 0]==4); 

end; 
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